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Abstract 

We have obtained a general solution of evolution equations for QCD twist-2 string 
operators in form of expansion over complete set of orthogonal eigenfunctions of 
evolution kernels in coordinate-space representation. In the leading logarithmic 
approximation the eigenfunctions can be determined using constraints imposed by 
conformal symmetry. Explicit formulae for the LO scale-dependence of quark and 
gluon twist-2 string operators are given. 
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Introduction 



Recently, there has been a renewed interest in QCD evolution of skewed parton dis- 
tributions [| . Skewed parton distributions play a crucial role in description of hard, 
exclusive QCD processes |T], ^|, £|, [| which are actively considered P, |7J, g, |9| as tools 
for investigation of new aspects of non-perturbative QCD dynamics [H| . However, it 
is clear that before a non-perturbative information can be reliably extracted from experi- 
mental data, all perturbative aspects, such as QCD evolution, have to be understood. So 
far, the main effort has been devoted to studies of evolution equations for skewed parton 
distributions in the momentum representation [[I], 0, |, [T^, [13|, [14], [15], [16], [IT], [TjJ. This is 
certainly the most natural choice in the forward, deep inelastic scattering limit, but, as ob- 
served recently [19|, as soon as there is a longitudinal momentum transfer between initial 
and final hadron states, the corresponding amplitude can be as conveniently represented 
in terms of momentum- as coordinate-space skewed parton distributions. 

In a recent paper |19| we have considered evolution of charge-conjugation odd and even 
flavor non-singlet coordinate-space skewed quark distributions in the leading-logarithmic 
(LO) approximation [20]. The resulting solutions have form of a Neumann series expan- 
sions J2lj over a set of Bessel functions of the first kind J%/2+j, with j assuming odd, 
respectively even integer values. That the Bessel functions enter the game is not surpris- 
ing at all - it has been known for a long time pOf that they form a representation of the 
conformal group in the coordinate-space and, as this symmetry is preserved at the LO 
22, they simultaneously diagonalize corresponding evolution kernels. What is more 



interesting is that the resulting expansions are valid uniformly for all values of the asym- 
metry parameter, in contrast to the situation in the momentum-space. There, one has to 
distinguish between cases where momentum fraction is smaller, respectively larger than 
the asymmetry parameter. A compact formula, based on the expansion in terms of local, 
multiplicatively renormalizable operators, exists only in the former case [§, 13, 15. 16. [0| . 



Evolution equations for twist-2 string operators have been known for a long time 
Written in the coordinate-space representation, they have been considered so far 



20, 24 



as an intermediate step in derivation of evolution equations for skewed parton distributions 
in the momentum space. However, given the transparent form of the solution to the 



problem of scale-dependence of coordinate-space skewed parton distributions found in |T0 



it is natural to ask whether similar formulae can be obtained for evolution of twist-2 string 
operators themselves. As we demonstrate below, this problem has a positive answer. The 
main goal of the present paper is to discuss an explicit solution to the evolution equations 
for string operators in a form of expansion over orthogonal eigenf unctions of corresponding 
evolution kernels. 

As we consider expansion of a non-local operator in terms of an orthogonal set of 
functions, it is natural that coefficients of such an expansion are themselves non-local 
operators. Their scale dependence is governed by a set of evolution equations which 
in general result in a mixing between coefficients corresponding to different orthogonal 
eigenf unctions. However, this term is present only when the analysis goes beyond the 
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leading order (LO) approximation. 

The remaining presentation is organized as follows. First, we introduce a general so- 
lution to QCD evolution equations for twist-2 string operators in a form of expansion 
over a set of orthogonal eigenfunctions of evolution kernels, find expansion coefficients 
and discuss their scale-dependence. Next, we perform an explicit analysis of the evolution 
equations in the LO approximation. The crucial point here is that appropriate sets of 
orthogonal eigenfunctions can be determined using constraints imposed by the conformal 
symmetry, .i.e. without solving the eigenvalue problems. In this approximation, the ex- 
pansion coefficients are identical with conformal string operators introduced in As a 
result, we obtain explicit expressions for QCD evolution of quark and gluon string oper- 
ators. Finally, we summarize. 



Evolution equations for QCD string operators 

Let us start from a discussion of a general structure of QCD evolution equations for string 



operators. As in Ref . ||20|| , we introduce the following definitions of quark, respectively 
gluon twist-2 operators 



0(a,/3) 



_,at + (3 . „ 



a + (3 a — (3 
z, z 



a - (3 



a + (3 a — (3 
z, z 



rib 
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In the above formula z is a light-like vector, z 2 
path-ordered exponential: 



0. The square brackets denote the 



rb 

[az,bz] = V exp[—ig / z^A^ytz) dt\ 

J a 



(3) 



which ensures gauge-invariance of the above definitions. Note that a z describes the 
center of the string composed from quark fields and the gluon line between them while 
(3 z corresponds to its "length" , understood simply as the difference between coordinates 
of quark or gluon fields at the string ends. 

Charge conjugation odd and even quark string operators are obtained from 0(a, 0) by 
taking its components symmetric, respectively antisymmetric in f3. Denoting the former 
and the latter by + (a,/3), respectively 0~(a,(3) we define 



+ (a,P) = -(0(a l ,/3) + 0(a l ,-P)) , 
0-(a,P)= l -(0(a,i3)-0(a,-i3)) . 



(4) 



Similarly to local operators, bare string operators have UV-divergences and should 
be renormalized. In complete analogy with the former case one can define renormalized 
string operator through the relation: 



R (u) = Z(u,u/)®O r (u/). 



(5) 
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Here we have introduced a compact notation ® for the two-dimensional integral 

/oo rco 
da / dp 
-oo JO 

and uj = {a,/?}. Function Z(ll>,u/) corresponds to a renormalization constant in a local 
case. We suppose that in MS-scheme it has the following structure: 

Z(u, uj') =5( W - W ') + E^E a k s Z nk {u, uj') , (6) 



n>l 6 ' k>n 



where a s = and e = 2 - d/2, d being the space-time dimension. Coefficients Z n k 

can be calculated in perturbation theory. In the MS-like schemes they do not depend on 
the renormalization scale \x. 

The renormalization group equation can be derived by taking a total derivative over 
fi of the left- and right-hand sides of (|5|). In this way one arrives at an evolution equation 
for the renormalized string operator Or(uj) |20], 24 1 

H-^0 R (u) = a s V(uj,u}')®0 R (u}'). (7) 

The evolution kernel V(u,u') is defined as 

a s V(u, uj') = h-^Z(oj, uj") ® Z~ X (J' , uj') . 

dfj, (8) 

Z~ x {uj" , uj') ® Z(u>', uj) = 5{uj" — uj) . 

Although it has not been indicated explicitly, V(uj,uj') is a function of the QCD coupling 
constant: 

vw) = E<(/^ (n W), (9) 

n>0 

except of the leading order (LO) approximation, when it is given simply by a function 
V^(uj,uj'). 

After these preliminaries, we shall construct now a general solution of the operator 
evolution equation flj]). To this end, let us assume that one can find a solution of the 
following system of equations for eigenfunctions and corresponding eigenvalues of the 
evolution kernel: 

V(uj, uj') ® tpi(oj') = 7#>i(u/) 

<Pi(u) (g> V(uj,uj') = ji<fi(uj) 

In the theory of integral equations the second equation is called conjugated to the first 
one [|25| . 



i.e. 



Eigenfunctions <Pi(ui) and <fj(ui) are orthogonal with respect to the scalar product ®. 

{<Pi, (fj) = uj>i(uj) ® (fj = SijNi , (11) 
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where iVj is some normalization constant. Note that beyond the leading order eigenfunc- 
tions <Pi(u), <fj(w) and eigenvalues ji depend on the running QCD coupling because the 
kernel V(uj,uj') itself becomes a function of a s (/x 2 ), see eq.(|9]). 
Solution of eq.(0) can be written in the following form: 

Or(u) = £ S„ Vn(w) • (12) 
n>0 

The operator- valued coefficients S n can be obtained by taking the scalar products of both 
sides of the above equation with eigenfunctions (p n \ 

S n = ^-<p n (J)®0 R {J) (13) 

As it follows, S n themselves are non-local operators. Moreover. eq.(0) leads to the fol- 
lowing set of evolution equations for operators S n : 

d d 

H~rS n + V S m {<f n , tp m ) = a s J n S n , (f> n = f3(a s )- ip n . (14) 

dfi m > da s 

The second term on the left-hand side induces mixing between operators corresponding 
to different orthogonal eigenfunctions. It is absent in the LO because eigenfunctions (p n 
do not depend on a s in this approximation. As a consequence operators S n become mul- 
tiplicatively renormalizable. This case will be discussed in more details in the next section. 

LO evolution of string operators 

In the following we will follow notation of |2D1 . To the leading-order accuracy, the evolution 
kernel is given and is therefore independent on a s (/x). Arguments 

based on the conformal symmetry can be used to find explicit solutions of the eigenvalue 
equation (|i~0p. In Ref . [p| the nonlocal conformal operators were introduced for the first 



time and the conformal symmetry arguments were used to find the solution of (0) in 
terms of an integral representation over a complex conformal spin j. Our goal here is to 



represent a solution in the form of an expansion (1121) ■ Is is equivalent to the former one, 
but much simpler from the numerical point of view. 

Now, let us consider solutions to the LO evolution equations in the framework set 
up in the previous section. First, note that the two-dimensional eigenvalue problem ([TO) 
can be reduced to a one- dimensional one by using translational invariance to separate the 
dependence on the variable a, which describes the position of the center-of-mass of the 
string. This is done by substituting 

^(a,/9) = e^hm, fa (a, (3) = e^/,(^) . (15) 

In this way, <fij(a,f3) and 4>j(a, 0) can be interpreted as wave functions of bound systems 
of two particles, with a momentum k and an internal motion described by fj(k/3) and 
fj(k/3), respectively. Let us first consider the case of flavor non-singlet quark operators. 
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Inserting ([15]) into the eigenvalue equations fllTf ) one arrives at the following system of 
equations p0| : 

rpi 

H lfj = / V{p 1 ,p 2 )f j {pi) = 7jfj(p2), 

J " _ _ (16) 

V(pi,p 2 )fj(p2) =lifj{Pl)- 



pi 



where p = k(3. The reduced kernel V(pi,p 2 ) is given by 



V(pi,p 2 ) = -35(pi - P2) + 



2 sin(pi - p 2 ) 4pi cos(pi - p 2 ) 



P2 



P2(pl - P2) 



+ A5{ Pl - P2 ) C j ■ (17) 



Equations ( ]T6| ) have the same form as the ERBL-equation for evolution of the pion distri- 
bution amplitude [22[ The latter is defined as the reduced matrix element of twist-2 
non-local quark operator between the pion state and the vacuum. Conformal symmetry 
was shown to play a crucial role in finding solution of the LO ERBL equation a long-time 



ago [22], p3| . In complete analogy, we shall now demonstrate how the symmetry con- 
siderations allow to determine eigenfunctions fj{kf3) and fj{k(3). Note that in the present 
case the evolution is considered at the operator level, without any reference to matrix 
elements. 

J 2 , the Casimir operator of the collinear conformal group SO(2,l), can be represented 
as an differential operator acting on the parameters a, (3. After substitution (|U|) one 
obtains |2( 



J 2 /,(p) ^i/,(p) 
J 2 f 3 (p)=L 2 f)(p) 



2 ^- 2 +^ + p 2 k 2 + 2 
dp 1 



dp 



fi(p) 



d 2 

2 — + n 2 k 2 



Up) 



Note that the explicit form of J 2 depends on whether it acts on fi(kj3) or on fi(k(3). 
Using equations (|T6| ) and ([18]) it is easy to check by an explicit calculation that J 2 indeed 
commutes with the 'Hamiltonians' Hi 2 , as required by the symmetry arguments. As it 
follows, the eigenfunctions of Hi 2 have to be the same as the eigenfunctions of L 12 : 

LJAP) = J(J + L 2 f 3 {p) = j(j + l)] 3 {p) , (19) 

where j are positive integers. As Li and L 2 are related by 

L2P 2 = P 2 U (20) 

one obtains the following relation between fj(p) and fj(p) 

fAp) = p 2 fj(p)- (21) 

Solutions of the set of differential equations (|19D are given by Bessel functions: 

f s {p) = Cj (k)p^ 2 Z 1+1/2 (p), Up) = -c 3 {k)p l ' 2 Z ]+ i /2 {p) , (22) 
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where Cj(k), Cj are yet undetermined normalization constants. As the operator Or(u) 
has a well-defined local limit as (3 — > 0, the validity of expansion ( |12|) requires that fj((3) 
are regular in this limit as well. In terms of fj and fj the scalar product ( [LTD reduces to 



(MM) = / dpf 3 (p)f 3 (p) = NAj , (23) 

<J 

and therefore fj and fj must be sufficiently regular such that the integral exists. These 
conditions are satisfied when we identify Zj + i/ 2 (p) with Jj + y 2 {p), the Bessel functions of 
the first kind. In particular, the scalar product ( p3|) reduces to the well known integral 
({t, j} = 1,2,3,...) [0: 

(fi(p)fj(p)) = c(k)iCj(k) f ^J 3+1/2 (p)J l+1/2 (p) = cik^ik)^ . (24) 

For reasons which will shortly become clear we set Cj(k) = k" 1 ! 2 and Cj(k) = k 3 ^ 2 . With 
this choice the nonlocal operators Sj defined according to (|13|) assume the form: 

/OO /"OO i 

dae lka / d(3J(3J 1/2+J (\k\P)0(a,(3)^. (25) 
-oo JO 

As noted above, they are multiplicatively renormalizable f20fl : 

^S(l/2 + j, k; p 2 ) = a slj S(l/2 + j, k; p 2 ), S(l/2 + j, k; p 2 ) = LjS(l/2 + j, k; p 2 ), 



7(i) 

L i = (a™)"" > 7(J) = iQQ — Cf (3 + TTJTfj — 4(*0' + 1) + 7E 



(26) 

Here, is the corresponding anomalous dimension known, e.g., from analysis of deep- 
inelastic scattering process. 

Now we are in position to write down the decomposition ([12]) explicitly: 

r°° dk 

3=1 



r°° dk ■? 00 

0(a,(3), 2 = / -e-^/3-^(l + 2j)^Ji/ 2+ ,(|fc|/5)5(l/2 + j,A;;^) (27) 

J—oo ATT -. 



It remains only to check the normalization. As p 2 = /1q and therefore Lj = 1 for all j we 
should have _ _ 

/OO POO 
da' dp"Jl(a,p\a',p')0(a,(3') (28) 
-00 Jo 

with 

/JjL OO 
—e^ a '-^ (3-1(3'^ J2(l+2j)J l/2+j (\k\(3)J 1/2+j (\k\(3') = 5(a-a')6((3-(3') 
i=i 

(29) 

To this end, recall that the series 

00 

prlfi £(1 + 2j)J 1/2+j (\k\(3)J 1/2+j (\k\(3') = W - (3') (30) 
i=i 
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does not depend on k and represents the Neumann expansion of a 5-function in terms of 
Bessel functions [19|, 21]. Equation (p9|) follows immediately. Note that this observation 



explains also the choice of the coefficients q, q. 

Equation (p7|) describes solution of the LO evolution equation as an expansion in terms 
of orthogonal eigenfunctions of the evolution kernel. On the other hand, it can be inter- 
preted as an expansion of a quark string operator in terms of conformal string operators. 
Series expansion in terms of Bessel functions is known in mathematical literature as the 
Neumann series, see [[H] for more detailed information. Evolution of flavor non-singlet, 
charge-conjugation odd and even quark string operators can be obtained from (|27f) by 
taking its symmetric, respectively antisymmetric in (3 components. 

Conformal operators (^5[) are defined only for positive, integer values of the conformal 
spin j. On the other hand, the solution found in |H| is based on an integral representation 



which requires analytical continuation of j into the complex plane. In our previous paper 
19|j we have explicitly checked that integration over j indeed reproduces the series ([27]) . 

Let us now turn our attention to operators with flavor singlet and positive charge 
parity quantum numbers. Here we have to consider mixing between gluon and quark string 
operators, see equations (|2|) and (^) for corresponding definitions. Note that the quark 
operator Q = 0~ receives contribution from all Np active flavors. An analysis analogous 
to the flavor non-singlet case leads to the following basis of the nonlocal conformal quark 
and gluon operators: 

/OO /"OO 
dae* ka / dP[3 3 / 2 J 1/2+j (\k\f3)G(a,l3)^, 
-OO JO ,„..s 

/oo roo i — \ J 

dae* ka / d(3pJ l/2+] {\k\(5)Q{a,l3)^. 
-oo JO 

Here j = 2 + 2n, n = 0, 1, 2, 3, ... is the conformal spin. Operators Sg(1/2 + j, k; jj?) and 
Sq(1/2 + j, k; n 2 ) depend on the renormalization scale according to 



S G {l/2 + j, k; /i 2 ) = [A+L+ - A_L_] S G {l/2 + j, k; fx 2 )+ 
+1gq {1 ^- [L+ - L_] S Q (l/2 + j, k; ti) , 



(32) 



S Q {\/2 + j, k; v 2 ) = [A+L_ - A+L_] S Q (l/2 + j, k; fi 2 ) + 
+-^- -i= [L + - L_] S G (l/2 + j, k; ^ 2 ) . 



(33) 



Here we have introduced a shorthand notation: 

7± 0) 

pog( P 2 /A 2 ) V~ 11. m 

7± = \ {iGG + Iqq ± ^~D) , D = (7 GG - 7 QQ ) 2 + 47 GQ 7 QG . 



(34) 
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In convention adopted in this paper the anomalous dimensions ^ggiIgqiIqg differ by 
a sign from those given in Ref.|20|. Note that although not explicitly indicated, j- 
dependence of the anomalous dimensions is implied here. 

To find the evolution of the string operators G(a, /3), Q(a, P) we have to rewrite them 
as series expansions in terms of the conformal string operators (|5T]) . Consider first the 
gluon operator G(a,/3). Using the Neumann expansion of the 5-function: 

oo 

PS(P-P') = Y;^v + 2 + An)J v+1+2n (\k\{3)J u+1+2n (\k\(3'), (35) 

n=0 

one easily obtains the desired formula: 

/OO /'OO 
da' / dp' 5 (a - a')5(p - P')G{a', P'\ li 2 ) = 
-oo JO 

/f dk 00 
J Z7Y n=0 



/ e- lto r 5/2 E( 5 + ^n)Js 2n (\k\f3)S G (l + 2n, k; fa 2 
J-oo Zn n=Q 



(36) 

Taking into account (|32]) one finally arrives at the following expression for the evolution 
of the gluon operator: 



r°° dk 00 
G(a,P;iu 2 ) = / -e^P^ 2 £(5 + An) Jt +2a (\k\0)x 



x { [X + L + - X-Ll j=2n+2) S G (I + 2n, k; /x 2 ) + (37) 



+1GQ- 



J(j=2n+2) W <-A 2 

2n+l 



D 



i L + ~ L -}(j=2n+2) + 2n ' fc 5 



Evolution of the quark-singlet operator can be found in a similar way. First, we 
establish an expansion in terms of non-local conformal operators 



r°° dk 00 

Q(a,/3;/i 2 ) = y ^_e-^r 3/2 E(5 + 4n)J |+2n (|fc|/3)5 Q (| + 2n,A;; yU 2 ). (38) 



dk 

n=0 

The resulting solution to the evolution equations reads 

<-°° dk 



/OO Wlp ' — ' 
-ika p-3/2 ( 5 + An) J 5 2n ( | k\P) X 

x |[A + L_ - X_L+} /j=2n+2) S Q (I + 2n, k; /i 2 ) + (39) 
+ (2^Tl) 7n [L + " L -h=2n+2) + 2n, fc; /i 2 )| . 
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Equations fl27|), fl37|) and (^) are new and represent the main results of this paper. 
As advocated in the previous section, they describe the scale dependence of QCD string 
operators in form of expansion over a set of eigenfunctions of corresponding evolution 
kernels. Evaluating matrix elements between appropriate nucleon states one finds the 
explicit form of scale dependence of skewed parton distributions in the coordinate-space 
representation Jl9|| . Due to translational invariance, the integral over k can be trivially 
performed and one is left with an expansion of corresponding coordinate-space skewed 
parton distributions in terms of orthogonal set of Bessel functions, valid for all values 
of the asymmetry parameter. Hence, the solution in the coordinate-space is much more 
transparent than in the momentum representation. 



Summary 

We have argued that solutions of evolution equations for twist-2 QCD string operators 
can be written as series expansions in terms of the orthogonal sets of eigenfunctions of 
evolution kernels. We have obtained a general form of evolution equations for operator- 
valued expansion coefficients and found that in a general case evolution results in a mixing 
between coefficients corresponding to different eigenfunctions. 

In the LO approximation the eigenfunctions of evolution kernels can be found explic- 
itly using conformal symmetry arguments and, at the same time, evolution equations for 
expansion coefficients become much simpler than in the most general case. As a result the 
expansion coefficients can be identified with conformal string operators. Those with flavor 
non-singlet or charge-parity odd quantum numbers are simply multiplicatively renormal- 
ized and those corresponding to flavor singlet and positive charge-parity are subject only 
to the usual mixing between gluon and quark operators. In both cases explicit formulae 
for scale dependence of quark and gluon string operators have been derived. 

Taking corresponding matrix elements one obtains immediately solution to the prob- 
lem of scale dependence of skewed parton distributions in the coordinate-space. Here, 
contrary to the situation in momentum space, the resulting expansion over a set of or- 
thogonal eigenfunction is valid for all values of the asymmetry parameter. 
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